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A microscopic model for the FeAs plane in the newly discovered Iron-based superconductors 
is proposed with the emphasis on the role of the orbital degeneracy between the Fe Zd^z and 
3dyz orbital in the crystal field environment. The model predicts commensurate antiferromagnetic 
ordered ground state for the parent compounds and a d-wave superconducting ground state at 
finite doping. Inter-orbital couplings plays an interesting role in enhancing both orders. Correlated 
variational wave functions are proposed for both ordered phases. 

PACS numbers: 



The newly discovered Iron-based superconductors F- 
doped LaOFeAspj] and its derivatives [3, |3|,y,|5|] with rare 
earth substitution arose many interests in the commu- 
nity. Similar to the high-T^ cuprates, the Iron-based su- 
perconductors also have a layered structure in which the 
role of the Copper oxygen plane in the high-Tc cuprates 
is played by the FeAs plane. In both systems, the charge 
dopant and the doped carrier are separated spatially 
which produces an ideal environment for the coherent 
motion of the latter at an incommensurate band filling. 
This similarity may imply that a still higher Tc is possi- 
ble in this series of compounds by optimize the interlayer 
distance and planeness of the FeAs layer. 

Band structure as calculated from Density functional 
theory(DFT) and its DMFT improvement had been ex- 
amined by several groups [Sl [Zl, la, Ifll, llO|. Not to one's 
surprise, most of the spectral weight close to the Fermi 
energy are contributed by the Fe 3d orbital, as the Fe is 
the only ion in the system that posses a unclosed shell. 
The DFT calculation also predicted a antiferromagnetic 
ordered ground state for the parent compound LaOFcAs 
with a ordered moment about 2.3 Bohr magneton, which 
may also be consistent with some preliminary experiment 
results on this system. 

Based on the results of the DFT calculations, several 
phenomenological theories for the doped Iron-based su- 
perconductors were proposed [ill [13, ll3|. In this paper, 
we follow instead the quantum chemistry considerations 
and propose a tight binding microscopic model for the 
FeAs plane. In our model, the orbital degeneracy be- 
tween the Fe 3^x2 and Zdyz orbital plays an essential 
role. Our model predicts commensurate antiferromag- 
netic ordered ground state for the parent compound and 
a d-wave superconducting ground state at finite doping. 
Inter-orbital couplings plays an interesting role in en- 
hancing both kind of orders. We also propose correlated 
variational wave functions for both phases. 

A key issue for the construction of a microscopic model 
for the FeAs plane is to elucidate the role of the orbital 
degeneracy of the Fe 3d orbital. In the presence of crys- 
tal field, the five fold degenerate Fe 3d atomic orbital will 



split according to the irreducible representations of the 
crystal symmetry and not all five 3d orbital are equally 
important for the low energy physics. According to the 
DFT calculation, the crystal splitting of the Fe 3d orbital 
is small and all five Fe 3d orbital should be included in 
the model in principle. However, the inclusion of the elec- 
tron correlation effect can enhance the crystal splitting 
and vice versa. Here we assume the crystal field splitting 
to be sufficiently large for the discussion of the low en- 
ergy physics and base our discussion solely on symmetry 
considerations. 

In the FeAs plane, each Fe ion sits at the center of in- 
version of a squashed(along the normal of the FeAs plane) 
tetrahedron formed by four neighboring As ions(see fig- 
ure 1). Using the lattice constants reported in the 
literature [7|, one estimate that the Fe-As bond make a 
angle of about 58.8 degree with the normal of the FeAs 
plane, while in a perfect tetrahedron the corresponding 
angle should be 54.7 degree. In such a local environ- 
ment of S4 symmetry [9|, the five fold degenerate 3d or- 
bital split into three nondegenerate(3(i3z2_r2, 3^3,2 _j,2, 
and idxy) orbital and one doubly degenerate orbital(3da;z 
and 2>dyz){see figure 2). If we assume that the crystal 
field is contributed mainly by the four neighboring As 
ions, one would expect the 3^3x2 _j.2 orbital to have the 
lowest energy, since the As tetrahedron in squashed in 
z-direction and the electron cloud of the Zd^z'^-j.^ orbital 
has the best chance to avoid the As ion(see figure 1). 
The next lowest energy orbital would be 2tdxy whose lobe 
points toward the vacancy between the As ions in the 
plane. If the tetrahedron is not squashed, the 'idxy or- 
bital would be degenerate with the Sd^^^-r^ orbital. The 
squashing distortion of the tetrahedron also lift the de- 
generacy between the remaining three orbital, namely 
Sdx2_y2,Sdxzi'idyz. Since the tetrahedron around the Fe 
ion is squashed toward the plane, the 3^3.2 _j,2 orbital 
should be higher in energy than 3dxz and Sdyz orbital. 
In the tetragonal structure of the FeAs plane, the degen- 
eracy between the 3dxz and 3dyz orbital is protected by 
symmetry. 

The divalent Fe ions in the FeAs plane have six elec- 




FIG. 1: The squashed tetrahedron formed by four neighboring 
As ion around each Fe ion in the FeAs plane. From the lattice 
constant reported in the Uterature, a=4.020 A, b=2.35 A, the 
angle (j> between the Fe-As bond with the normal of the Fe- 
As plane is estimated to be 58.8 degree, while in a perfect 
tetrahedron the angle would be 54.7 degree. Here we use a 
coordinate system in which the x and y axis point along the 
As-As bond directions. 



trons in their 3d shell. According to the above splitting 
scheme, one would expect a very peculiar situation to 
occur in the parent compounds in which the Fermi en- 
ergy lies within two degenerate bands which are both half 
filled. We think this is the reason that the parent com- 
pounds differ from usual band insulator or band metal. 
It may also hold the key for the lightly doped LnOFeAs 
system to show such remarkable properties as having a 
Tc as high as 50 K^. 

In the following, we assume the crystal splitting to be 
large enough to neglect the filled 3^322 .^2 and "id^y band 
and the empty 3d^2_y2 band in the discussion of low en- 
ergy physics. With this simplification, we are left with 
a model with two degenerate bands. In this model, each 
Fe site can accommodate at most four electrons in the 
two degenerate bands. In the parent compound, each Fe 
site have two electrons on average in the two degenerate 
bands and the system is thus half filled. The interactions 
between the electrons on the same Fe site can be classified 
into three types, namely the intra-orbital Hubbard re- 
pulsion U , the spin independent part of the inter-orbital 
Hubbard repulsion Ui, and the spin dependent part of 
the inter-orbital Hubbard repulsion J (the usual Hund's 
rule coupling). From the definition of these terms, one 
easily see that that the inequality U > Ui > j > 
should be satisfied. These interaction terms is repre- 
sented by the model Hamiltonian 



+C/i 2j("»,2:z,T + ni^xz.,i){ni,yz,] + ni,yz,i) 

i 

i 

in which ni^xz,i — cJa;^ tCi.^z,! denotes the number 
of up spin electron in the Sd^^ orbital. Si^^z = 
^a,f3Ci^xz,f3 denotes the spin density on the Sd^z 



2 t,xz,a 

orbital. The meaning of other terms in the equation is 
similar. 




FIG. 2: The five d-orbital in the crystal field of the tetra- 
hedron formed by the four neighboring As ions around each 
Fe ion. The 3d^2_y2 and the Sd^y orbital are shown here 
in a projection onto the xy plane. In the projection graph, 
red(green) filled circles denote As ion above(below) the xy 
plane. Here we only show the 3dxz orbital, the 3dyz orbital 
can be obtained from the 3dxz orbital through a rotation of 
■| along the z-ax;is. 



In the FeAs plane, the Fe ions are too remote from each 
other to have direct hoping and the As ion is needed to 
bridge the hoping between the Fe ions. The As ion has 
three 3p orbital in its outmost shell. Since the 3pz orbital 
is far away from the Fermi energy, we consider only the 
hopping path mediated by the 3px and 3py orbital of As 
ion. Symmetry considerations shows that the hoping me- 
diated by the As 3px and Spy orbital conserve the orbital 
index, namely, an electron initially on the Fe Sdxzi'idyz) 
orbital can only hop onto the 3dxz{3dyz) orbital of desti- 
nation Fe ion(the inclusion of the hopping path mediated 
by 3pz orbital will break such a symmetry) (see figure 
3). Another peculiarity about the As ion aided effective 
hoping is that the hoping integral between next-nearest- 
neighboring Fe ions is anisotropic (the hoping between 
nearest neighboring Fe ions remains isotropic). This 
anisotropy depends on the orbital and is a reflection of 
the anisotropic nature of the 3dxz and Sdyz orbital. Tak- 
ing all these considerations into account, we arrive at the 
following model for the effective hoping between the Fe 



Ht = -t 






'^i+S,xz,a + '^i,yz,a'^i+5,yzMl + -W.C. 
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o-''i+5" ,2;z,cr "•" ^i,yz,(T'^i+&' ,yz,a) 



,(2) 



in which S = e'x ot e^ denotes the nearest neighboring 
vectors on the square lattice form by Fe ions. S — e'^ + e'y 
and 5 = e^ — e^ denotes the next nearest neighboring 
vectors. Since there are two hoping paths between the 
nearest neighboring Fe ions but only one for the next- 
nearest neighboring Fe ions, t is roughly a factor of two 
larger than both ii and t2- From symmetry consider- 
ations, one can also show that all of the three hoping 
integrals, t, ti and ^2 can be taken positive with a suit- 
able choice of gauge. The dispersion relation of the 3dxz 
band is given by 



Ex 



-2t{cos{kx) + cos{kx)) — 2ii cos(fc2; 



-2t2Cos{kx - ky), (3) 



the dispersion relation of the 3dyz band can be obtained 
from that of 3dzz band by the reflection in momentum 

The total Hamiltonian of the system is the sum of 
the on-site term Hu and the hopping term Hf. Before 
looking at the phase diagram of the model, we first es- 
timate the model parameters. According to DFT calcu- 
lations, the Hubbard repulsion U and Ui are about AeV 
in magnitude and the Hund's rule coupling J is about 
leV in magnitude. For the hopping terms, as both the 
Fe^dxz — AsSpa, separation and the hoping integral be- 
tween the two are of the order of leV, we estimate the 
effective hoping integral between the Fe ions mediated 
by the As ions to be also of the order of leV. Thus the 
system has a moderate level of electron correlation and 
moderate level of frustration(as indicated by the ratio 
between t and ti). 

Now we discuss the possible phase diagram of the 
model. In the half filled parent compounds, each site is 
occupied on average by two electrons. Since the intra- 
orbital Hubbard repulsion U is larger than the inter- 
orbital Hubbard repulsion f/i, the two electrons tend to 
occupy different orbital and the remaining Hund's rule 
coupling tends to align up the spin of the electrons on 
both orbital. Thus each individual Fe ion carries ap- 
proximately a unit spin angular momentum and shows 
two Bohr magneton of magnetic moment. However, the 
inclusion of the hoping terms between the Fe ions will 
induce antiferromagnetic spin exchange between neigh- 
boring Fe ions within each band. Since the frustration 
is moderate, one can expect an commensurate antiferro- 
magnetic order to appear in the ground state of parent 




FIG. 3: A schematic representation of the As 3p orbital aided 
lioping between neighboring Fe ions on the square lattice. 
Note that the As ions lies above or below the Fe plane. Sim- 
ple symmetry consideration shows that the electron on Fe 
3dxz{3dyz) orbital can hop only through the As 3px{3py) or- 
bital. This explains the conservation of orbital index and 
the anisotropy of the hoping Hamiltonian. The electron wave 
function on the yellow(blue) lobes have positive(negative) val- 
ues. 



compound. As a result of the Hund's rule coupling, the 
antiferromagnetic order in the Sd^z and the Sdyz band 
enhance with each other. With this understanding in 
mind, we can write down the following variational wave 
function for the ground state of the half filled system. 



I*) = |SDW),, (g, |SDW)i 



(4) 



in which ISDW)^^ and |SDW)j;;, denote the SDW ordered 
state in the Sd^z and the 3dyz band. 

In our model, the two bands are correlated with each 
other through the inter-orbital Hubbard term and the 
Hund's rule coupling. To describe such correlation effect, 
we introduce the following Jastrow-type variational wave 
function 



I*>=5V5j|SDW),,0|SDW), 



(5) 



in which < g,gi,g,j < 1 are Gutzwiller factors in- 
troduced to describe the correlation effect caused by 
intra- and inter-orbital Hubbard resulpsion and Hund's 
rule coupling. The operators D, V and S are given by 
D = Hi{ni,xz,]ni^xz,i + "-i,y2,T"»,yz,i)' ^ = I]i('^i,x2,T + 
ni^xz,i){ni,yz,] +nuyz,i) and S = Y.i{Si,xz •S'i,yz + |). In 
the antiferromagnetic ordered ground state, we expect 
the correction due to the correlation effect to be small as 
the main effect of the inter-band correlation is embodied 



in the mutual enhancement of the AF order in the two 
bands. 

At finite doping, the antiferromagnetic exchange 
within each band wiU induce pairing instabihty in the 
d-wave channel as in the high-Tc problem. Unlike the 
high-Tc problem, now we have two degenerate bands 
which are correlated with each other by the inter-orbital 
Hubbard repulsion and the Hund's rule coupling. As in 
the antiferromagnetic ordered state, it is likely that the 
inter-orbital coupling will induce mutual enhancement of 
the pairing instability in both bands. However, the inter- 
orbital coupling will inevitably increase the effective mass 
of the Cooper pair in each band and thus reduce the su- 
perfluid density. With this understanding in mind, we 
can write down a correlated variational wave function 
for the superconducting state as we have done for the 
antiferromagnetic ordered state. 

I*) = .Af 5,/|d - BCS),, |d - BCS)^„ (6) 

in which |d — BCS)xz and |d — BCS)j,2 denote the d-wave 
BCS superconducting state in the 3dxz and the 3dyz 
band. We note that the Gutzwiller factor describing the 
Hund's rule coupling, g'J, now plays a very similar role 
as that of the projection operator in the construction 
of the matrix product type ground state of the AKLT 
modelil4]. Thus the superconducting state can be more 
appropriately called a condensate of entangled Cooper 
pairs rather than independent Cooper pairs. 

The arguments presented above are rather crude and 
a more careful study of the model proposed in this pa- 
per is obviously needed to elucidate its relationship with 
the superconductivity observed in F-doped LnOFeAs. In 
principle, the arguments can be made more quantitative 
through systematic mean field treatment or Gutzwiller 
approximation. A direct numerical calculation on the 
variational wave function is also possible. However, be- 
fore any such effort is made, it is important first to know 
the actual value of the crystal field splitting by, for ex- 
ample, the optical measurement, and check the validity 
of neglecting the Sd^z^-r^, Sdj:2_y2 and the 3dxy orbital 
in the low energy subspace. 

In summary, a microscopic model for the FeAs plane of 
the newly discovered Iron-based superconductors is pro- 
posed. In our model, the orbital degeneracy between the 
Fe 5dxz and the Sdyz orbital plays an important role and 
may explain the peculiarities of the parent compounds 
of this family of new superconductors. The model pre- 



dicts a commensurate antiferromagnetic ordered ground 
state for the parent compound and a d-wave pairing su- 
perconducting ground state in the doped case. In the 
latter case, the inter-orbital coupling transform the sys- 
tem from a condensate of independent Cooper pairs into 
a condensate of entangled Coopers and it would be in- 
teresting to see effect of inter-orbital couplings on the 
dynamical properties of the new superconductor. 
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